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An analytical result for renormalization of the jump of the heat capacity AC/Cn by the anisotropy 
of the order parameter is derived in the framework of the model proposed by Haas and Maki [Phys. 
Rev. B 65, 020502(R) (2001)], for both prolate and oblate anisotropy. The graph of AC/Cn versus 
the ratio of the gaps on the equator and the pole, Ae/Ap, of the Fermi surface allows a direct 
determination of the gap anisotropy parameter Ae/Ap using data from specific heat measurements. 

PACS numbers: 74.20.Rp, 74.25. Bt, 74.70.Ad 


In a recent Rapid Communication, Haas and MakiB 
discuss a model for the gap anisotropy in MgB2 , a mate- 
rial which has attracted a lot of attentfon from condensed 
matter physicists in the past two years. A central issue 
in this work is to propose an analytic model for analyz- 
ing the thermodynamic behavior. Assuming a spherical 
Fermi surface a simple gap anisotropy function is sug- 
gested, A(k) oc 1 -I- az^, where z = cosd, and 6 is the 
polar angle. This angular dependence is a natural linear 
combination of spherical harmonics (Legendre polynomi- 
als Poiz) and P2{z)), and can therefore be interpreted in 
terms of an s+d model. 
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FIG. 1: Jump of the specific heat AC/Cn as function of ratio 
between "equatorial" and "polar" order parameter Ae/Ap. 
For one and the same jump AC/Cn we have prolate and 
oblate solution. 

This model leads to useful results for the temperature 
dependence of the upper critical field Hc2 and of the spe- 


cific heat, which can be fitted to the experimental data, 
thereby determining the optimal anisotropy parameter 
a. Note that a — Ap/A^ — 1, with Ap = A(z = 1) and 
Ae — A{z = 0). For example, the dimensionless ratio 
AC/Cn of the specific heat jump at Tc relative tolhe 
normal state specific heat, measured bypWang et al.p is 
a tool for estimating a. Haas and MakitI focused on the 
case a > 0. 

For oblate-shaped gap anisotropy (Ag > Ap), in- 
dicated b}j__recent experiments on the critical field 
anisotropy,cla we take a < 0. The sign of a is relevant 
for this physical property, but irrelevant for the jump 
of the heat capacity. Note that for a < —1 the model 
gives a gapless superconductor, which is not applicable 
for MgB2. Following the weak-coupling BCS approach 
by the authorsEJ we derived the analytic expression 


AC 

Cn 


12 l + 4a/3 + 38aV45 + 4a3/15 + 0^/25 


7C(3) 
1 

Ae/Ap 


4a/3 + 6a2/5 

Ae/Ap = 


-4a3/7- 
1 

1 + a 


ay 9 


(1) 


valid for arbitrary a. This expression can be used in 
experimental data analysis and may lead to two solutions 
(oblate, a < 0; prolate, a > 0), for a given specific heat 
jump. The relevant example is shown in Fig. |l|. 

For the jump AC/Cn ~ 1-18 considered by the au- 
thors,EJ we recall a = 1 and contribute a = —0.6877. The 
corresponding anisotropics are Ae/Ap = 1/2 (prolate) 
and Ae/Ap = 3.202 > 1 (oblate solution). 
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